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Abstract. 

We here present results of three-dimensional Smoothed Particle hydro and magnetohydrodynamics simulations of rotating jets, 
also including the effects of radiative cooling, precession and velocity variability. Using initial conditions and parameters which 
are particularly suitable for the DG Tau microjet, we have been able to approximately reproduce its complex knotty morphology 
and kinematics. We have also obtained radial velocity maps which are in good agreement with the data obtained by Bacciotti et 
al., thus indicating that their interpretation that the DG Tau microjet is rotating is correct. Finally, we have found that a magnetic 
field of the order of ss 0.5 mG is sufficient to collimate the jet against the lateral expansion that is caused by the centrifugal 
forces. 
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1. Introduction 

The seek for the real nature of the accretion-ejection mecha- 
nism underlying the formation of the jets is one of the major 
contributions of the detailed analysis of Herbig-Haro (HH) jet 
data. The statement that MHD models can be responsible for 
generating the jets associated with Young Stellar Objects is still 
a matter of debate. The theoretical discussion of the so-called 
magneto-centrifugal models for jet launching (e.g., Blandford 
& Pavne [T982l Konigl 1982, Spruit [T996l . can be strongly con- 
strained by the existence of observational data. 

Recent observations of microjets associated with T-Tauri 
stars, namely, DG Tau, RW Aur, TH28 and LkHa321 
(Bacciotti et al 120021 Coffey et al l2004t . have revealed impor- 
tant trends in the jet radial velocity field that are consistent with 
a rotation pattern inside the jet beam. The presence of rotation 
could be understood in the context of the magneto-centrifugal 
models as due to jet launching from a Keplerian accretion disk 
along with the magnetic field lines which are anchored into the 
disk-star system and provide the collimation of the beam. 

Motivated by these recent findings, we have carried out a 
set of fully three-dimensional simulations of rotating jets, also 
taking into account the presence of radiative cooling, jet pre- 
cession, and velocity variability (to allow the production of in- 
ternal knots). All these are necessary ingredients to address the 



Send offprint requests to: A.H. Cerqueira 



kinematical complexity of the microjets and we are here par- 
ticularly interested in the DG Tau microjet case, which is one 
of the best studied in the literature. We have also investigated, 
in few models, the collimating effects of magnetic fields upon 
rotating jets. 

2. Numerical method and the simulated models 

In order to examine the behaviour of the microjets and their 
early evolution near the source under the effects of both ro- 
tation and precession, we here employ modified versions of 
our three-dimensional Smoothed Particle Hydrodynamics (3-D 
SPH) and Magnetohydrodynamics (3-D SPMHD) codes which 
were previously designed to investigate the large scale struc- 
ture and evolution of protostellar jets either in the absence (see, 
e.g., de Gouveia Dal Pino 2001) or in the presence of magnetic 
fields (e.g, Cerqueira & de Gouveia Dal Pino 2001). 

Our computational domain is a rectangular box that mim- 
ics the ambient medium and has dimensions -25Rj < x < 25Rj, 
and -lORj < y,z < lORj, where Rj is the initial jet radius (which 
is also the code distance unit). The Cartesian coordinate sys- 
tem has its origin at the center of the box and the jet flow is 
continuously injected into the bottom of the box [at the inlet 
which has coordinates r = (— 257?j, 0, 0)] with top-hat density 
and pressure profiles. The jet and gas pressures are assumed to 
be in equilibrium at the inlet. Inside the box, the SPH particles 
are initially regularly distributed in a cubic lattice. An outflow 
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Table 1. The jet models 



Model 


Tpu\ 


^"prec 


rotation? 


B-field? 


A 




8 years 


no 


no 


B 


8 years 


8 years 


no 


no 


C 


8 years 


8 years 


yes 


no 


D 






yes 


no 


E 






yes 


yes 




Fig. 1. Midplane velocity field distribution for the jet models A 
(top), B (middle) and C (bottom), at tju = 2.5. See the text for 
details. 



boundary condition is assumed for the boundaries of the box. 
The particles are smoothed out by a spherically symmetric ker- 
nel function of width h, and the initial values of h were chosen 
to be OARj and 0.2Rj for the ambient and jet particles, respec- 
tively. The adiabatic index of the ambient medium and the jet 
is assumed to be y = 5/3, and an ideal equation of state is 
used. The radiative cooling, due to collisional excitation and re- 
combination, is implicitly calculated using a time-independent 
cooling function for a gas of cosmic abundances cooling from 
T 10 6 to 9 x 10 3 K. 

The adopted parameters for the simulations are appropri- 
ate to the conditions generally found in protostellar jets and, in 
particular, in the DG Tau microjet. We take a number density 
ratio between the jet and the ambient medium 77 = «j/n a = 10; 
nj = 2000 crrT 3 ; an average ambient Mach number M d = 
Vj/c a = 19, where Vj - 300 km s is the average jet velocity of 
DG Tau (e.g., Bacciotti et al. 2002, Pyo, Kobayashi, Hayashi et 
al. 2003), and c a 16 km s _1 is the ambient sound speed for a 



gas with an average temperature of 10 4 K; and Rj = 5.6 x 10 14 
cm (see e.g., Dougados et al. 120001 Bacciotti et al. 2002). 

We have carried out several simulations considering both 
continuous and pulsing jets (see below). In the latter cases, we 
have adopted a sinusoidal profile to describe the ejection veloc- 
ity time- variability at the inlet: v (t) = Vj[l + A ■ sin(^f)L 
where A is the velocity amplitude and r pu i is the period of 
the oscillation. The values for these parameters were obtained 
from observations of the DG Tau microjet (see, e.g., Lavalley 
et al. 1997, Raga et al. UMT\ which suggest A = 0.33, and 
T" P ui = 8 years = 0.71?^ (where % = ^j/c a ~ 11.3 years cor- 
responds to the transverse jet dynamical time). In the numer- 
ical models where the jet precession has been also taken into 
account, we have assumed an equal precession period (i.e., 
Tprec = Tpui = 8 years) and a precession half-angle = 5° 
as indicated by the DG Tau jet observations (Lavalley-Fouquet 
et al. 2000, Dougados et al. l2000l 

For the computation of the jet rotation around its main axis, 
we have assumed that the flow conserves its angular momen- 
tum after its passage through the Alven surface (i.e., v$-r = con- 
stant, where is the toroidal velocity in the normal direction to 
the jet axis). This is consistent with magneto-centrifugal mech- 
anisms for jet launching (see, e.g., Spruit 1996 1 and is also sug- 
gested by the observations of the DG Tau microjet (Bacciotti et 
al. 2002 1 which imply = (4.362 x 10 2<) jr) cm s _1 where r is 
the radial distance to the jet axis in cm. This gives v<s - 8 km 
s _1 at the jet surface (r = Rj = 5.6 x 10 14 cm), and « 55 km 
s _1 near the jet axis (at r - 0.15/?j) It should be emphasized 
that we are here mostly interested in carrying out simulations 
of jet regions near the source but beyond the Alfven surface 
(e.g., Spruit 1996 1 where the flow has accelerated to velocities 
larger than the Alfven speed and is believed to be already colli- 
mated, and which are, at least in some cases, partially resolved 
by optical observations (as for instance, DG Tau). 

Table 1 summarizes the properties of the simulated models. 
In model A, we have run a precessing jet with constant velocity 
at injection and no rotation. In model B, we have considered a 
pulsing precessing jet also without rotation, and in model C we 
have run a jet which is simultaneously pulsating, precessing, 
and rotating. All these three models are purely hydrodynam- 
ical and share the same initial conditions as described above. 
Their velocity distribution is compared in Figure 1 after they 
have propagated over a time f/fj = 2.5. In model C, the jet 
is precessing counter clockwise while rotating clockwise. The 
jet of model A develops a V-shape structure at the head which 
is commonly detected in precessing jets (e.g., Masciadri et al. 
2002] Masciadri & Raga l2"0TJ3l Rosen & Smith l2"U0ll and is 
due to the impact of the parcels propagating in different direc- 
tions in the wiggling flow with the bow-shock front. In the jet 
of model B, the combined effects of pulsation and precession 
produce internal working surfaces propagating in different di- 
rections. In this case, a narrower bowshock structure develops 
at the head that makes it to propagate slightly faster than the 

1 We notice that the code does not allow, by construction, any SPH 
particle to occupy the jet axis, so that the equation above for has 
no singularities, as it is employed only for values of r which are larger 
than zero. 
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Fig. 2. (a) Midplane density contour (left panel) and velocity 
field (rigth panel) distributions for the model C (see Table 1), 
at t/t& = 1 which approximately corresponds to the dynamical 
age of the DG Tau microjet; (b) Average radial velocities across 
the jet for the different positions along the beam that are marked 
in Figure 2a, above. We notice that the jet axis is actually lo- 
cated at z = O.lRj due to the effects of the precession. The 
vertical bars indicate the maximum and minimum values of the 
radial velocities obtained for each radial distance (the distances 
are in units of the jet radius); (c) Radial velocity differences be- 
tween positions across the jet which are symmetric with respect 
to the jet center, as a function of the radial distance. Crosses: 
our results. Full line: an approximate fitting to Bacciotti et al. 
(2002) observations and Pesenti, Dougados, Cabrit et al (2004i 
results (the distances are in units of AU and the velocities in 
units of km s ). 

jet of model A. The rotating jet of model C essentially bears 
the same features of the previous models but develops an even 
larger head than model A due to the centrifugal forces that tend 
to push the jet flow outwards in the radial direction (see also 
Figure 3 and discussion below). 

Let us now concentrate on the early evolution of the rotat- 
ing, pulsing, precessing jet of Figure lc (model C) whose char- 
acteristics qualitatively resemble those of the DG Tau microjet 
(Dougados et al. 2000 Bacciotti et al. 2002 1. Figure 2a shows 
the jet of Figure lc at an earlier time that corresponds to the 
inferred dynamical age of the DG Tau microjet (~ 1 1.3 yr; see 
e.g., Dougados et al l20001 Lavalley et al l 19971 Pyo, Kobayashi, 
Hayashi et al. 2003). The density contour plot presents nearly 
similar features to those observed in the DG Tau microjet: an 




Fig. 3. Midplane velocity field distributions for the purely hy- 
drodynamical model D (left) and for the MHD model E (right) 
of Table 1, at t/td = 1. The distances are in units of Rj 38 
AU). 

elongated, highly collimated beam near the jet base, a slight 
sideway displacement of the beam that has been previously in- 
terpreted as an evidence for precession of the jet axis, and a 
prominent bow shock (at « -8Rj in Figure 2a) which can be 
compared with the Bl knot of DG Tau (Dougados et al. 2000 
Lavalley-Fouquet et al. 2000 1. 

For a comparison with the radial velocities obtained from 
the HST observations of DG Tau by Bacciotti et al. (2002|l, in 
Figure 2b we have built diagrams of the radial velocity across 
the jet (taken at an angle of 45° with respect to the jet axis, in 
order to mimic the DG Tau inclination with respect to the line 
of sight) in four different slices along the flow which are sep- 
arated by 0.5Rj - 20 AU and are marked in Figure 2a. [These 
correspond to the same slices examined in Figure 1 of Bacciotti 
et al.: x = -24.5^ (Position I), -24Rj (Position II), -23.5R } 
(Position III) and -23.0Rj (Position IV).] Figure 2c depicts, for 
each diagram of Figure 2b, the radial velocity differences be- 
tween positions which are symmetric with respect to the central 
region of the jet, as a function of the distance from the jet axis. 
We notice that the differences in radial velocity are negative, 
therefore indicating that the net velocity is oriented towards the 
—y direction. These negative shifts in velocity are obviously 
a consequence of the fact that the jet of Model C is rotating 
clockwise. As expected, these shifts are null in the jets that are 
only precessing, like those of Figures la and b. The resulting 
radial velocity shifts in Figure 2c, Av r < -15 km s , are very 
similar to those obtained by Bacciotti et al. which are approxi- 
mately reproduced in the figure by the solid lines. We note that 
these solid lines also reproduce the results recently obtained by 
Pesenti, Dougados, Cabrit et al. (2004i assuming that the DG 
Tau microjet is launched magnetocentrifugally by a warm ac- 
cretion disk. The consistency between the present results and 
their predictions is an additional support for the magnetocen- 
trifugal mechanism for jet launching. 

Finally, aiming to check the collimation effects of the mag- 
netic fields in rotating jets, we have performed simulations of 
continuous, non-precessing flows under the presence of heli- 
cal fields. Given the present uncertainties related to the real 
orientation and strength of the magnetic fields in protostellar 
jets we have assumed, as in previous work (e.g., Cerqueira & 
de Gouveia Dal Pino 2001) a force-free helical magnetic field 
which also extends to the ambient medium and whose radial 
functional dependences for the toroidal and longitudinal com- 
ponents are respectively given by: 
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W = B [f(rW' 2 , (la) 

B x (r) = B [l-f(r)(r + d)] 1 ' 2 , (lb) 

where r = \Jy 2 + z 2 is the radial distance from the jet axis, 
f(r) = Q.5Cr 2 /(r + Q.5d) 3 , and the constants C and d are as- 
sumed to be 100 and 3/?j, respectively. In these equations, Bo 
gives the maximum strength of the magnetic field and corre- 
sponds to the magnitude of the longitudinal component at the 
jet axis and the magnitude of the toroidal component at r — 3Rj 
(see Cerqueira & de Gouveia Dal Pino 2001). 

As an example, Figure 3 compares the early evolution of 
one of these MHD runs (model E of Table 1) with a purely 
hydrodynamical jet (model D of Table 1). Both have the same 
initial conditions of the previous models which are appropriate 
to DG Tau microjet. The MHD model has an initial value of the 
ratio between the magnetic and the gas pressure B = /? ma g/Ah = 
B 2 /(8np) =10 near the jet surface which corresponds to a mag- 
netic field intensity ~ 0.5 mG. The centrifugal force upon the 
fluid in the radial directions causes the beam radius of the hy- 
drodynamical jet to increase as the central parts of the flow tend 
to move laterally outwards. Though less obvious, we have pre- 
viously noticed the same unbalancing of radial forces and jet 
enlargement in the rotating jet of Figure lc (Model C), partic- 
ularly in its head. When a helical magnetic field (with surface 
P =i 10) is introduced in the system, we find that the tension 
force (or hoop stress) associated to the toroidal component of 
the magnetic field is able to collimate the jet. Smaller values of 
the surface ft were found to be unable to prevent the lateral ex- 
pansion and similar results were also obtained for pulsing jets 
under the same conditions (see discussion below). 

3. Discussion and conclusions 

In the present work, we have explored the role of rotation on 
the evolution of YSO microjets addressing, in particular, the 
DG Tau microjet that has recently shown strong observational 
evidence for rotation (Bacciotti et al. 2002 1. With the help of 
3-D hydro and magnetohydrodynamical simulations, we have 
investigated several models that also included the effects of 
precession and pulsation of the flow. Our main purpose here 
was to verify through fully 3-D simulations whether the sig- 
nature of rotation was really unique and unmistakable in the 
jet flow when other effects, like precession and pulsation, were 
also considered. The main conclusions of this study could be 
summarized as follows: 

1. The morphology and kinematics of the DG Tau microjet 
was best reproduced in our simulations by model C (Figure lc 
and 2a), which presents a jet with a sinusoidal velocity variabil- 
ity (with mean velocity of 300 km s _1 , half-amplitude of 100 
km s _1 , and a period of 8 yr), a precession of the outflow axis 
(with half-angle of 5° and 8 yr period), and a rotation around 
its axis with a velocity 8 km s _1 at the jet surface, and 
as 55 km s _1 near the jet axis,. These latter conditions sat- 
isfy angular momentum conservation and are consistent with 
magneto-centrifugal jet launching models. As remarked before, 
these characteristics approximately match those suggested by 
the observations reported in the literature. We note also that 



previous numerical studies of the morphology and kinematics 
of this object (e.g., Raga et al. !2001i have obtained similar con- 
clusions although they have not taken into account the effects 
of rotation. 

2. The simulated radial velocities for model C have revealed 
good agreement with the data obtained by Bacciotti et al 2002 
therefore suggesting that their interpretation that the DG Tau 
microjet is rotating must be correct. 

3. The inclusion in the models of a helical magnetic field 
permeating the jet and the ambient medium, with a maximum 
value of p =i 10 near the jet surface (corresponding to a toroidal 
magnetic field ~ 0.5 mG), has provided a collimation of the 
outflow which otherwise would tend to expand laterally while 
propagating downstream (Figure 3). Though the MHD mod- 
els here investigated were only simple cases involving rotat- 
ing continuous (or pulsing) jets without precession, the results 
above suggest that the presence of magnetic fields of the order 
of ~ 0.5 mG in the DG Tau microjet would be sufficient to col- 
limate the beam against the centrifugal forces. In fact, a simple 
estimate of the magnetic tension forces that would be required 
to balance the centrifugal forces at the DG Tau microjet sur- 
face gives a magnetic field strength B = (Anpv 2 ^) 1 ' 2 ~ 0.2 mG 
(where p n/nyi ~ 3.3 10 21 g ctrT 3 and = 8 km s _1 ), which 
is compatible with the required value above from the simula- 
tions. We further notice that this would be the expected strength 
of the magnetic field at radial distances ~ 40 AU from the 
source. Magnetic flux conservation would then imply strengths 
~ 100 G near the stellar surface which are compatible both 
with the values expected for protostellar sources and also with 
magneto-centrifugal mechanisms for jet production. 
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